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We investigate the electronic transport properties of unbiased and biased bilayer graphene nanorib- 
bon in n-p and n-n junctions subject to a perpendicular magnetic field. Using the non-equilibrium 
Green's function method and the Landauer-Biittiker formalism, the conductance is studied for the 
cases of clean, on-site, and edge disordered bilayer graphene. We show that the lowest Hall plateau 
remains unchanged in the presence of disorder, whereas asymmetry destroys both the plateaus and 
. conductance quantization. In addition, we show that disorder induces an enhancement of the con- 

D ' ductance in the n-p region in the presence of magnetic fields. Finally, we show that the equilibration 

of quantum Hall edge states between distinctively doped regions causes Hall plateaus to appear in 
[■ — , the regime of complete mode mixing. 

i— I ' 

PACS numbers: 73.63.-b, 72.80.Ng, 81.05.ue, 05.60.Gg 

^ ' I. INTRODUCTION 

In recent years, after the success of fabrication of both monolayer— and multilayer graphene sheets 2 - - — , there has been 
a lot of interest in the transport properties of graphene nano-ribbons, especially the behavior of low-energy charge 
carrier excitations^—. Ideal monolayer graphene is a gapless semimetal with zero density of states at the Dirac points. 
The low-energy electronic excitations in the vicinity of the Dirac points have linear dispersions, and are described 
by an effective massless Dirac Hamiltonian. The low-energy electrons in bilayer graphene, on the other hand, have 
a quadratic dispersion relation. Both for monolayer and bilayer graphene, the wave functions are composed of two 
' sublattices A and B, and give rise to the chirality of the charge carriers. Therefore, the charge carriers in graphene 
are chiral and has awaken an enormous interest in graphene, for instance with regards of effects such as the Berry 
O ■ phase, which is 7r in monolayer graphene and 2-k in bilayer graphene. Although intrinsic bilayer graphene is a zero-gap 
O \ semimetal, it exhibits very interesting properties when a gate voltage is applied, which makes bilayer graphene into a 
tunable band gap semiconductor^. The band gap determines the threshold voltage and the on-off ratio of field effect 
y—i ' transistors and diodes. Therefore, bilayer graphene is more convenient for applications in nano-electronic industry 
^ , than monolayer graphene^. 

One of the exotic phenomena that has been observed in monolayer and bilayer graphene is the anomalous quantum 
Hall effect — The nature of massless chiral Dirac charge carriers in monolayer graphene gives rise to this property 
of the Hall plateaus, that behave as a xy = ±<tq(N + 1/2), with N being the Landau level index and uo = 4e 2 //i. 
The factor of 4 originates from the valley and spin degeneracies. In undoped bilayer graphene, the sequence of Hall 
plateaus, with a xy = ±<7qN, were observed. The first plateau at N = is missing which implies that bilayer graphene 
is metallic at the neutrality point, while the standard quantum Hall effect in bilayer graphene can be recovered by 
applying a gate voltage. The quantum Hall states, fully quantized due to the presence of a magnetic field, as well 
as broken-symmetry states at intermediate filling factors such 0, ±1, ±2 and ±3, were experimentally observed by 
Fcldman et al^. 

In a perfect nanoribbon, the electron transmission via subbands due to lateral confinement of the electronic states 
implies the quantization of the conductance in units of Go = 2e 2 /hJ£ Recently, the zero-temperature conductance 
d • of free-disordered monolayer and unbiased bilayer graphene nanoribbons in the presence of a uniform perpendicular 
magnetic field was calculated^. The conductance in monolayer graphene nanoribbon is given by 2(n + l/2)Go for 
the case of zigzag edges, and uGq for the case of armchair edges. On the other hand, it was shown that in a bilayer 
graphene nanoribbon the conductance is quantized as 2(n + l)Go for zigzag edges, and nGo for armchair edges, where 
n is an integer. 

The quantum Hall effect and quantized transport in graphene junctions in the bipolar (p-n), and unipolar (n-n 
or p-p) regimes was investigated theoretically and experimentally by several group s 15 ' 16 . Long, et al^ by using the 
Landaucr-Buttikcr formalism, showed that on-site disorder induces the enhancement of the transport in monolayer 
graphene p-n junctions in the presence of a magnetic field. On the other hand, they showed that in the n-n junction, 
the lowest plateau survives in a sufficiently broad range of on-site disorder strengths. They also showed that in a 
particular range of disorder strengths new plateaus {i.e. G = 3e 2 /h and e 2 /h) emerg o 15 ' 17 , something also observed 
experimentally. 



in 



o 



- 1—1 

X 

5-H 



FIG. 1: (Color online) (Left): Schematic 3D view of bilayer graphene incorporating all coupling energies. (Right): Schematic 
picture of zBGNR with length M and width N atoms, and the definition of the principal layer. 



Transport measurements in high quality bilayer graphene pnp junctions have also been performed, and electron 
mobilities up to 10000 cm 2 /(V s) have been measured for gapless systems, and an on-off ratio up to 20000 for gapped 
systems^ Moreover, the fractional-valued quantum Hall plateaus due to equilibration of quantum Hall edge states 
between distinctively doping regions have been observed Consequently, the conductance exhibits plateaus arising 
from the mixing of edge states at the interfaces. 

In this paper, we investigate the conductance of the chiral massive carrier in the presence of a uniform perpendicular 
magnetic field, both for unbiased and biased bilayer graphene nanoribbons configured as n-n and n-p junction a 18 ' 19 . 
In addition, we study the influence of on-site and edge disorders on the conductance. For these purposes, we use the 
tight-binding model and Landauer-Biittiker formalism together with the non-equilibrium Green's function approach^. 
It should be noted that, in general, the absolute magnitude of the magnetic filed can cause reflection at the boundaries 
of the electronic devices^. We restrict our attention in this article to the case of a system in which the reflection on 
the boundaries due to the magnetic field can be ignored. 

The paper is organized as follows. In Sec. II, we introduce our model and formalism, such as the tight-binding 
Hamiltonian for bilayer graphene junction, and a recursive method for calculating the Green's function. In Sec. Ill, 
our numerical results for the conductance of disordered and biased bilayer graphene junctions in the presence of 
magnetic field are presented. Finally, we conclude in Sec. IV with a brief summary. 



We consider a bilayer graphene ribbon with Bernal stacking (AB)22, as a conductor connected to the left and the 
right leads as depicted in Fig. Q] The model Hamiltonian is given by, 



where H center , Hl and Hr are the Hamiltonian of the center region, the left and the right leads, respectively. Two 
leads arc considered as perfect semi-infinite bilayer graphene nanoribbons. Notice that the leads are also structured by 
the Bernal stacking bilayer graphene. We consider the nearest-neighbor tight-binding Hamiltonian with one 7r-orbital 
per site on the lattice. The effective one-body Hamiltonian of bilayer graphene in the presence of the magnetic field 
is given as follows, 



II. MODEL AND METHOD 
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where a\ i and ai.i (bj i and are the creation and annihilation operators at sublatticc A (B) in the layer 1=1,2 
at the zth site, respectively. The intralayer nearest-neighbor hopping energy is 70 = 3.16eV, the hopping energy 
between on-top sublattices A and B in different layers is 71 = 0.39eV and furthermore, 73 = 0.315eT^ denotes the 
hopping energy between not on-top sublattices A and B between two layers^. Another hopping energy between the 
nearest-neighboring layers, 74 = 0.04eV , is very small compare to 70 and can be ignored. In the presence of the 
external perpendicular magnetic field B, the hopping integral acquires the Peierls phase factor given by e l ^ iJ where 
<pij = f? A ■ dl/0o, with the magnetic flux quantum 0o = h/e. We use the Landau gauge as A = (—By, 0,0). The 
applied magnetic field is considered to be only on the center region. Vij reduces to the bias voltage El (E r ) on the 
left (right) lead and can be controlled by the gate voltage. The electrostatic potential changes from the right lead to 
the left lead and is assumed to be linear as vi = k(E R — El)/(M + 1) + El, k = 1,2, M, where M is the length of 
the center region (see Fig. 1). We consider on-site disorder Wi, being a random variable with a uniform distribution in 
an interval [-W/2, W/2] with the disorder strength W which exists only in the center region. The size of the central 
region, i.e. conductor, is given by AN x M atoms. Here we define the asymmetric between two layers, A, indicating 
the difference between on-site energies. The current can be calculated from the Landauer-Buttiker formula^ as 



2e 

T 



deT LR (e)[f L (e) - f R (e)} (3) 



where / Q (e) = l/(cxp[(e — eV a ) /ksT] + 1), (a = L,R), is the Fermi distribution function in the leads. To calculate 
the transmission coefficient Tlr{(), we use 

T LR (e) = Tr\r L GT R &] (4) 
where the line width function r Q , that describes the coupling between the conductor and the leads is given by 

r a (e) = i[£;(e) - £»(e)], (5) 
in which the retarded Green's functions is written as 

G ( e ) = u \^~T\ v^TT ( 6 ) 

C — ^center ~ ^L\ e ) ~ ^R^J 

The retarded self-energy ££,, due to the coupling to a-th lead is 

X r L = h LC (e-H L )- 1 h LC (7) 
Z r R = h CR (e - H R )-^ CR 

where h,Lc(hc R ) is the hopping Hamiltonian from the left lead to the center region (from center region to the right 
lead) and G a = (e — iJ Q ) _1 can be calculated by an iterative method numericall y 23 ' 24 . We assume an infinite 
stack of principal layers with the nearest-neighbor interactions. A principle layer is defined as the smallest group 
of neighbouring atoms planes such way that only nearest-neighbour interactions exist between principle layers ( see 
Fig. 1). Thus, we can transform the original system into a linear chain of the principal layers. By using this approach, 
we write the matrix elements of (e — H)G = 1 in the following form 2 ^ 

(e— H o)Go,o = 1 + HoiGi : o (8) 
(£ — ^00)01,0 = H Q1 Go,o + Hq\G2,o 

(e — Hoo)G n .o = H^Gn-ifl + HqiG u +i,q 

in which Hqo and Hqi describe the coupling within the principal layer and the adjacent principal layers, respectively. 
For simplicity, we assume that -ffoo = Hn = H22--- and Hai = i?i2 = H23-... Notice that G nm is the matrix element 
of the Green's function between the principal layers. It is easy to obtain an iterative set of equations for G raj o 

Gfi.O = iiG n _2 i ,0 + tiGn+2 i ,0 (9) 



for n > 2*. where 



k = (1 - U-iU-i - U-iU-i)' 1 ^-! (10) 
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where 

to = (e-Hao)- 1 !!^ (11) 
to = (e — #00) 1 H Q i 

then we can write 

Gi.o = t Goo + t G2 t o (12) 
= (to + toti)Goo + tiG4fi 

= (to + toil + ••■ + to...t n -it n )Go,o + t n G 2 i+i.o 

We solve Eq. ((9]) iteratively. This process is repeated until t„ + i,t„ + i < e. in which e is a tiny value and is chosen as 
small as one pleases. Therefore G 2 "+i j o — 0, and the transfer matrix is thus given by T = to + toti + toti^2 + ••• + 
to...t n -\t n . Accordingly we can write G\p = TGo,o and Go.o = TGifi. The self-energies of the conductor-leads are 
Sl = H^f and Y, R = H m T. Finally, the zero-temperature conductance G = limy_>o can be obtained by using 
the Landauer-Biittiker formalism as2£ 

G=^-T LR (e F ). (13) 
h 



III. NUMERICAL RESULTS 



In this Section, we present our numerical results for zero-temperature conductance of unbiased and biased zBGNR 
in the presence of a magnetic field as well as various types of disorders. We assume that the width of the nanoribbons 
is N = 45, the voltage of the left lead is El = —0.2. All the energies are in units of 70. We neglect the effect of 
the Zccman splitting and the spin-orbit interaction, which are important only at very low energies, at which disorder 
effects normally dominate.— We define dimensionless magnetic field as <f> = (3\/3 / 4)a 2 B j (f>a , where 2<f> is the magnetic 
flux in a honeycomb lattice. 

In connection with the formation of the Hall plateaus, it is necessary to consider ribbons with width greater than 
the magnetic length scale, Ib = yhfeB. We consider = 0.01 corresponds to Ib ~ 15A which is smaller than the 
considered ribbon size with width L y (N = 45) w lOnm, and length L X (M = 21) ps 5.5nm. 

The conductance of unbiased clean zBGNR as a function of Er is shown in Fig. 2 for various sizes, both in the 
absence (<j> = 0) and presence (<j> = 0.01) of a magnetic field. In the absence of a magnetic field, the conductance in the 
n-n region (En < 0) is quantized due to transverse confinement of the ribbon, and well described by G = 2(n + l)Go, 
with the minimum conductance of a zBGNR being 2Gq^. Moreover, the conductance is independent of the ribbon 
length at low Er values. As one can see in the n-n region for = 0, the energy spacings between plateaus are not 
equidistant (whereas in monolayer graphene they are) , because of the quadratic dispersion relation. The widths of 
the conductance steps are related to the energy scale between the successive modes in the energy spectrum. Therefore 
the conductance is sensitive to El values, and the number of plateaus increases with increasing bias voltage, i.e. 
\El — Er\. Also for Er < El there are no plateaus, identical to the case of monolayer graphene^. In the n-p 
region, Er > 0, the conductance occurs due to the chiral charge carriers tunneling between n and p regions, and the 
conductance is always less than the corresponding plateau value in the n-n region. In this region the conductance 
decreases with increasing length of the ribbon M, since the number of scattering centers increases. 

The effect of a high magnetic field, (f> = 0.01 on the conductance of a clean graphene junction in unipolar and 
bipolar regimes is shown in Fig. 2. In the n-p region the length dependence of the Peierls phase factor gives rise to 
a non-monotonic behavior of the conductance as a function of the length and energy Er, noticeably at very low Er 
values. This behavior is in contrast with the result obtained for the zero magnetic field. Our results show that the 
conductance of the clean sample in bipolar regime is suppressed dramatically in the presence of the magnetic field. 
However, the first and the second Hall plateaus survive (i.e. G/Gq = 2 and 4) in the n-n region. 

We also study the effect of asymmetry between two layers, A, i.e. when the two layers have different on-site energies. 
The asymmetry here leads to the opening of a gap between the conduction and valance bands. In Fig. 3, we plot the 
conductance of a free-disordered zBGNR as a function of Er in the absence of the magnetic field. In the n-n region, 
asymmetry leads to a decrease of the conductance, while in the n-p region, asymmetry results in an enhancement of 
the conductance. This effect can be described based on the channels of the charge carriers. In the n-n region, electrons 
are only charge carriers while in the n-p region, because of the existence of asymmetry between two layers, one layer 
is n-doped and the other is p-doped. Accordingly, both electrons and holes play a role in the transport. In the n-n 
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region, the conductance fluctuations occur for low \Er\ and increase with length size M. Importantly, opening a gap 
affects the transversal confinement and the quantized steps are destroyed by asymmetry in the n-n region. 

The conductance of a zBGNR for various A is plotted in Fig. 4. The asymmetry leads to an increase of the 
conductance in the n-p region, and a decrease in the n-n region. The conductance fluctuations increase with increasing 
A values. In the right panel of Fig. 4 we show the effect of the magnetic field on the conductance in the presence of 
asymmetry. In the n-n region the Hall plateaus are destroyed by asymmetry and accordingly the conductance reduces. 
On the other hand, in the n-p region asymmetry leads to increasing conductance. 

We are now in the position to introduce some disorder and study the effects of disorder on the conductance of 
zBGNR in the presence of a magnetic field B. We consider the effect of on-site disorder on the conductance with a 
uniform disorder distribution in the range of [-W/2, W/2]. In Fig. 5, we plot the conductance as a function of E R 
for the various W for = (left panel). The number of realizations is 200. In the n-n region, the conductance is 
suppressed. In the n-p region at low En values, the conductance increases when the strength of disorder increases. 
The conductance is independent of disorder in the large positive Er regions. 

We also consider a disordered zBGNR in the presence of the high perpendicular magnetic field. The quantum Hall 
effect in gapless bilayer graphene occurs at the filling factors v = ±1, ±2, ±3..., in which there are \u\ = \noh/eB\ edge 
modes propagating in the opposite directions at v > and v < 0. Here uq is the charge density. The conductance 
plateaus for unipolar regime is given by G nn /Go = G pp /Gq = 2(min(|i/£,|, = 2, 4, 6... and for a bipolar system the 

conductance is G pn /GQ = ^jj^qq^j) = 1, f , §, 2... in the two-terminal ohmic regimei 16 ' 18 It is worthwhile mentioning 
that those expressions can also be applied to monolayer nanoribbons in the ohmic regime since these two materials 
have similar resistivities and thus also similar mean free paths^. 

For the unipolar regime the edge states which are common between the left and right regions propagate between 
the two leads, while \vl — vr\ states do not contribute to the conductance. In the bipolar regime, on the other 
hand, mode mixing occurs in the interface of two regions and for complete mode mixing the conductance plateaus 
obey the aforementioned formula for G pn /Go- In unipolar regime our numerical calculations show that (0 = 0.01, 
right panel) the lowest Hall plateaus remain unchanged in the presence of small disorder strengths, whereas the Hall 
plateaus are destroyed with increasing disorder strengths. We see perfect Hall plateaus in the unipolar regime with no 
equidistance in the scale of Er. In the clean bipolar regime the Hall edge states are separated for electrons and holes, 
and leads to suppression of the conductance. In addition, small strengths of on-site disorder induce the enhancement 
of the conductance of zBGNR in the presence of a magnetic field in the n-p region. At small strengths of disorder, 
W < 1, in the n-p region, the conductance is enhanced due to the mixture of electron and hole edge states. Thus 
in the bipolar regime, mode mixing at interface leads to two- terminal conductances. On the other hand, for large 
values of disorder strength, the system enters the insulating regime and the conductance is very small for all El and 
Er. We expect that the lowest Hall plateau survives only within certain range of disorder strengths. In the inset 
of Fig. 5 the conductance is shown in the ohmic regime and it obeys G pn /Go = 2( [j^Tn^~i ) = lj §, |,2..., the edge 
state equilibration condition. It should be noted that the ohmic behavior has been observed experimentally in bilayer 
graphene pnp junctions as well as graphene p-n junction a 16 ! 18 . Our numerical results are in excellent agreement with 
the recent experiment. 

We also investigate the effect of asymmetry in the presence of on-site disorder and our results are shown in Fig. 6. 
The Hall plateaus are destroyed by asymmetry even for very low disorder strengths. Also, asymmetry destroys finite 
size quantization of the conductance in the n-n region. The strong fluctuations vanish in the n-n region for small Er 
values and in the presence of disorder, as is shown in Fig. 6 for the case of = 0. 

Another type of disorder which is indispensable in real nanoribbons is edge disorder—. This type of disorder is 
generated by eliminating carbon atoms randomly along the edges of GNR. Note that because of our limitations in 
this approach, we consider only one layer depth edge disorder. In Fig. 7, we have shown the effect of edge disorder on 
the conductance of zBGNR as a function of Er for two different nanoribbon lengths M. The conductance increases 
in the n-p region at low Er as compared to the clean system, and it is independent of disorder in large Er. It is 
important to investigate the persistence of the Hall conductance plateaus versus edge disorder. As we have shown in 
Fig. 7, the Hall plateaus remain unchanged in the presence of edge disorder as well as on-site disorder for = 0.01. 

IV. SUMMARY AND CONCLUSION REMARKS 

In summary, we studied the effect of on-site and edge disorder on the conductance of biased zigzag bilayer graphene 
nanoribbon subject to a uniform perpendicular magnetic field. Our approach was based on the non-equilibrium 
Green's function method and Landaucr-Biittikcr formalism. Our results show that the lowest Hall plateaus can 
survive in the presence of a broad range of disorder strengths in the n-n region, while an asymmetry between two 
layers destroys them. On the other hand, disorder induces an enhancement of the conductance in the presence of the 



6 



magnetic field in the n-p region. In addition, the conductance is enhanced due to asymmetry in the n-p region. We 
also showed that the Hall plateaus appear due to equilibration of the quantum Hall edge states in the different regions 
with electron and hole type charge carriers. 

Our approach can be extended to long-range disorder due to charge impurities, and also to the case of spin 
dependence of the electronic transport with ferromagnetic-gate in bilayer graphene nanoribbon sheets. 
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FIG. 2: (Color online) Conductance of a clean zBGNR as a function of Er for various lengths with El = —0.2 for A = 0.0 at 
(j> = and <j> = 0.01. 
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FIG. 3: (Color online) Conductance of a clean zBGNR with the asymmetry between two layers, as a function of Er for various 
lengths with El = —0.2 at zero-magnetic field. 




FIG. 4: (Color online) Conductance of a zBGNR as a function of Er for M = 21, El = —0.2, and various A values at (j> = 
and <j> = 0.01. 
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FIG. 5: (Color online) Conductance of a zBGNR as a function of Er for M = 21, El = —0.2, and various on-site disorder 
strengths at cj> = and <j> = 0.01. 
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FIG. 6: (Color online) Conductance of a zBGNR as a function of Er for M — 21, El = —0.2, A = 0.1 and various on-site 
disorder strengths at <f> = and <j> = 0.01 . 




FIG. 7: (Color online) Edge disordered conductance of a zBGNR as a function of Er in comparison with the clean system in 
the various values of the lengths for El = —0.2 at = and cj> = 0.01. 



